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OBJECTIVES 


One of the fundamental difficulties in modeling cloud fields is the large 
variability of cloud optical properties (liquid water content, reflectance, 
emissivity). The stratocumulus and cirrus clouds, under special consideration 
for FIRE, exhibit spatial variability on scales of 1 km or less. While it is 
impractical to model individual cloud elements, the research direction is to 
model a statistical ensembles of cloud elements with mean-cloud properties 
specified. The major areas of this investigation are: 1) analysis of cloud field 
properties; 2) intercomparison of cloud radiative model results with satellite 
observations; 3) radiative parameterization of cloud fields; and 4) development 
of improved cloud classification algorithms. 

RESULTS OF COMPLETED WORK 

The following results have been completed under the support of this 
grant. These results have been published in refereed journals or presented at 
science conferences. 

1 ) Parameterization of Stratocumulus Cloud Fields . A radiative 
parameterization scheme of stratocumulus cloud fields is developed based on 
our Monte Carlo simulations. The parameterization includes Rayleigh scattering 
and atmospheric absorption above the cloud layer. Zero surface reflection is 
assumed below the cloud, roughly simulating an ocean or lake surface. The 
parameterization is valid for cloud cover ranging up to 100%, for zenith angles 
less than 72.5° and for cloud optical depths between x = 3 and x = 49. 

2) Structural Characteristics of Marine Stratocumulus Cloud . It is 
found that stratocumulus cloud size distributions obey a power law, that cell 
horizontal aspect ratio increases with cell diameter, that these clouds are bi- 
fractal in nature, and that they appear to be homogeneous over regions of 
about 100 km x 100 km. These are important constraints in our modeling of 
these cloud fields. 

3) Pattern Recognition Method for Cloud Classification . Texture has 
been found to be a very powerful descriptor of cloud field structure. A new 
texture-based pattern recognition method has been developed which is a hybrid 
of Max-Min statistics and Generalized Co-occurrence Matrices (GCM). Using 
these techniques, it is found that cloud fields may be classified to within 95% 
accuracy using high spatial resolution single channel near-infrared data. 

4} Effects of Broken Cloudiness Upon Surface Reflectance . Monte 
Carlo simulations have been performed as a function of cloud cover, cloud 
optical depth, solar zenith angle and surface albedo. The results show that 
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broken cumulus cloudiness of about 20% cloud cover may not be detectable 
over reflecting surfaces with albedo as low as 0.25 - 0.30. Over high albedo 
surfaces with albedos greater than about 0.40, partial cloudiness may actually 
decrease the system albedo. 

5) Comparison of Radiative Transfer Theory with Observations by 
Satellite and Aircraft . Observations of cirrus and altocumulus clouds during 
FIRE by Landsat and the King Air have been compared to theoretical models of 
cloud radiative properties. Good agreement is found between observations and 
theory when water droplets dominate. Poor agreement is found when ice 
particles dominate. 

6) Cumulus Cloud Spacing and Clustering . Detailed observations of 
cumulus cloud size distributions, inhomogeneities, nearest-neighbor relation- 
ships and cloud field scales of clustering have been determined for ten cloud 
fields using Landsat data. 

7) Cloud Base Height Determination . Cloud base heights are critical 
for accurate radiative energy balance studies. A Hough Transform approach is 
used to estimate cumulus cloud base height from hjgh spatial resolution 
Landsat data. The approach employs a variety of image processing techniques 
to match cloud edges with their corresponding shadow edges. Cloud base 
height then is estimated by computing the separation distance between the 
corresponding Generalized Hough Transform reference points. Another method 
for determining cloud base height of cirrus clouds also was developed. The 
cross correlation between a transparent cloud and its shadow is used to 
determine their separation distance. It is estimated that cloud base height 
accuracies of 50-70 m may be possible using HIRIS and ASTER instruments. 

8) Automated Detection of Jet Contrails . AVHRR split window 
imageries are used to examine cirrus and contrail signatures. A two step 
approach is used for this purpose. A preliminary algorithm subtracts the 

1 1 .8 mm image from the 10.8 mm image to enhance contrails. Then a three- 
stage algorithm searches the difference image for the nearly-straight line 
segments which characterize contrails. 

9) Cloud Classification Using Texture Analysis . We use textural 
characteristics of cloud and ice covered surfaces as a tool of cloud 
classification in polar regions. Neural network technique, an artificial 
intelligence classifier, has been used by us as an approach to cloud 
classification. We also examine the loss of cloud classification accuracy as a 
function of spatial resolution by degrading the imagery through progressive 
averaging. Significant improvement in cloud classification accuracy can be 
obtained using 1 /2 -km spatial resolution data rather than the current 1 km 
resolution data available today from AVHRR and GOES. Cirrus classification 
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accuracy is especially compromised as the spatial resolution is degraded. 
However, the use of texture measures defined at the combination of pixel 
separations d = 1 , 4 improves classification accuracies by several percent even 
for 1 km spatial resolution data. Cirrus accuracy is significantly improved by 
use of multiple distance features. 

1 0) Developing Monte Carlo Radiative Transfer Code Under IDL . The 
Monte Carlo radiative transfer code has been rewritten in IDL. The radiances 
calculated by the Monte Carlo radiative transfer code has been verified against 
the plane-parallel spherical harmonics and discrete ordinate models. The Monte 
Carlo radiances are obtained in discrete angular bins. Since not all bins have 
the same number of photons, the radiances agree better for bins with more 
photons (higher radiance) and poorer for bins with fewer photons (lower 
radiance). This behavior is characteristics of Monte Carlo simulations. A 
number of simulations have been run, most notably to check Kobayashi’s 
results. We have run the case of a plane parallel lower cloud layer with a 
checkerboard pattern on top. Our flux results are nearly identical to those 
reported by Kobahashi. 

1 1 ) Retrieval of Cloud Microphvsics . Using a radiative transfer model 
to interpret AVHRR radiances, we have developed a practical method to 
retrieve effective droplet radius for liquid water clouds on a global scale. We 
performed sensitivity tests to investigate the possibility and the accuracy of 
retrieving effective droplet radius and developed a method to remove the 
thermal emission contribution to Channel 3 radiances. A methodology has also 
been developed to assess instrument noise of all AVHRR channels under 
operational conditions. We have performed a global survey of the effective 
droplet radii of liquid water clouds. 

12) Inferring Cloud LWP and Validation of Cloud Droplet Retrievals . 

We developed a method using visible, infrared and near infrared bands to 
retrieve cloud LWP in liquid water clouds using the simple relationship between 
LWP, cloud optical thickness and effective droplet radius. Cloud optical 
thickness is obtained from visible band and effective radius is retrieved from 
visible, near-infrared and infrared radiances. Note that this analysis only 
retrieves effective radius from the top of clouds, so the results may 
overestimate or underestimate LWP if r 9 is not a good approximation of 
vertically mean effective radius over the whole cloud layer. 

13) Coreaistration of AVIRIS and TIMS Imagery and Retrieval of Cloud 
Properties . A technique was developed for coregistering AVIRIS and TIMS 
imagery using both line by line adjustments and 2D scene adjustments. It was 
found that the 3-band ratio method of Gao and Goetz (1991) applied to AVIRIS 
imagery was superior in detecting cloud pixels on the cloud sunside whereas 
the thermal IR bands of TIMS were superior in detecting pixels on the 
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antisunside. It was also shown that the shadows cast by the short-level and 
relatively fast moving FWC clouds induced an average 1 °C temperature 
difference. 

14) Cumulus Three-Dimensional Cellular Structure . LANDSATTM 
channel 6 images are used to obtain the statistics. A cell recognition algorithm 
is developed to automatically examine the cloud cellular structures. This allows 
both unicellular and multicellular analyses. A linear least squares method is 
employed to model the shape of cloud cells with quadric surfaces. Clouds 
smaller than 1 km in diameter generally are unicellular and have smaller fractal 
dimensions. The larger clouds are mostly multicellular and have larger fractal 
dimensions. The cloud cells exhibit characteristics similar to those of the small 
unicellular clouds, with the exceptions: (1) that the cell size distribution 
appears to be better modeled by an exponential distribution and (2) that the 
larger cells have smaller horizontal aspect ratios, meaning that they have 
increasing tendencies to be more circular. For the 3000+ cells examined, it is 
found that 80% can be best approximated by the shape of a hyperboloid of one 
sheet. About 15% of the cells are hyperboloids of two sheets with convex 
domes in the positive z direction. The average number of cells in a cloud is 
found to increase with increasing cloud size slightly faster than linearly. A 
cloud with effective diameter of 5 km is composed of approximately 5 cells. 
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